Band Structure Engineering of Multinary Chalcogenide Topological Insulators 
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Topological insulators (TIs) have been found in strained binary HgTe and ternary I-III-VI2 chal- 
copyrite compounds such as CuTlSe2 which have inverted band structures. However, the non-trivial 
band gaps of these existing binary and ternary TIs are limited to small values, usually around 10 
meV or less. In this work, we reveal that a large non-trivial band gap requires the material having a 
large negative crystal field splitting Acf at top of the valence band and a moderately large negative 
s — p band gap Eg~^. These parameters can be better tuned through chemical ordering in multi- 
nary compounds. Based on this understanding, we show that a series of quaternary I2-II-IV-VI4 
compounds, including Cu2HgPbSe4, Cu2CdPbSe4, Ag2lIgPbSe4 and Ag2CdPbTe4 are TIs, in which 
Ag2HgPbSe4 has the largest TI band gap of 47 meV because it combines the optimal values of Acf 
and E^-P. 

PACS numbers: 73. 20. At, 71.15.Dx, 71.18.-l-y, 73.61.Le 



The search for new topological insulators (TIs) has in- 
tensified recently due to their scientific importance as a 
novel quantum state and the associated technological ap- 
plications in spintronics and quantum computingP, Q. 
So far, experimental realizations have been limited to 
a few classes of simple materials, including zinc-blende 
based HgTe quantum wellsfl-Q, Bii_2;Sba; alloys[l, 01 
and binary tetradymite semiconductors such as Bi2Se3 
and Bi9Te?if8l-[l0t. Most recently, the search for TIs 
has extended to ternary com pounds dMI, e.g., strained 
Half-Heusler compounds, in the hope that the presence 
of more chemical elements would bring greater material 
flexibility. Despite the success of identifying these TIs, 
the design of new TI materials with the following advan- 
tages is still desired: (i) reahzing a topological insulating 
state with a significant non-trivial band gap (i.e., larger 
than kT at room temperature) at its natural equilibrium 
state (i.e., not under external strain), (ii) easy integration 
with electronic and spintronic devices based on tetrahe- 
dral semiconductors, and (iii) easy to be synthesized or 
already have been synthesized. 

Based on the direct evaluation of the Z2 topological 
invariant, Feng et al.[l^ proposed that a series of I- 
III- VI2 chalcopyrite compounds (such as CuTlSe2) could 
have topologically non-trivial band structure, and some 
of them can realize a topological insulating phase in 
their natural equilibrium structure. This is an impor- 
tant observation because the chalcopyrite structure is de- 
rived from the zinc-blende structure, and the band struc- 
ture properties are well understood, mostly for solar cell 
applications [H, Some of the proposed Cu and Ag 
based TIs, such as CuTlSe2 and AgTlTe2, have already 
been synthesized experimentally [l^ However, the 

predicted band gaps of these TIs are very small, usu- 
ally around 10 meV or less, similar to that observed in 



strained HgTe[4|. 

In this Letter, we show that the non-trivial band gaps 
of zinc-blende derived compounds with inverted band 
structure are mainly determined by the crystal field split- 
ting Ac p at the top of the valence band and the size of 
the inverted s-p band gap Eg~P, which can be better 
tuned by changing the component elements in a multi- 
nary ordered compounds. A large non-trivial band gap 
requires the material having a large negative Acf and 
a large negative E^'^p as long as it has no band crossing 
at the Fermi energy. For I-III-VI2 topological insula- 
tors, because the band inversion requires the group-Ill 
elements to be large and heavy, whereas a large nega- 
tive Acf requires group-Ill elements to be small and 
light, the possibilities for obtaining a large TI band gap 
are limited. Through further cation mutation 19], large 
negative A,;^^' and Eg~P is achievable in quaternary II2- 
II-IV-VI4 compounds. We have identified four topologi- 
cal insulators (Cu2HgPbSe4, Cu2CdPbSe4, Ag2HgPbSe4 
and Ag2CdPbTe4), in which Ag2HgPbSe4 has the largest 
TI band gap of 47 meV. In the following, we will discuss 
the evolution of the band structure of zinc-blende derived 
structures and explain what kind of band structure can 
lead to the largest TI band gap. 

For a normal zinc-blende semiconductors such as 
CdTe, the band gap is between the s-like conduction band 
minimum (CBM) Fg^ state and the p-like valence band 
maximum (VBM) Tgy state, as shown in Fig. 1. The non- 
trivial band structure of a TI is characterized by the band 
inversion in the Brillouin zone(6l fisj. i.e., the position of 
the conduction and valence bands is switched. In zinc- 
blende compounds, the band inversion means that the 
Fgc level falls below the Fgu level. In the inverted band 
structure, the Fgc level is occupied, while the quadruply- 
degenerate Fg^ level is half occupied, making the Fermi 
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FIG. 1: (Color online) (a) The conduction and valence band 
splitting of cubic and tetragonal semiconductors, (b) A plot 
showing how the band structure of normal semiconductors 
transfers into the inverted and topological insulator band 
structures. Note that the subscript v (c) represents the state 
belongs to the valence (conduction) band in the normal band 
structure. 



level stay at the Fg^ level and the system become a zero- 
gap semi-metal. This is the case for bulk HgTe. 

To open a band gap and change the zinc-blende semi- 
metal HgTe into a topological insulator, one has to induce 
a crystal field splitting Acf by reducing the Td symmetry 
of the zinc-blende structure to, e.g., I?2d, by applying an 
epitaxial strain or forming a quantum well|^]. For D2d 
symmetry, the half-filled Fg„ state splits into Fg^ and Fy^ 
state, and a gap can be opened around the occupied Fy^ 
(Fgi,) and unoccupied Fg,, (^Tv) levels (Fig. 1(b)) if Acf 
is positive (negative). On the other hand, the crystal 
field splitting can also be induced by chemical ordering, 
e.g., by mutating two Hg (group II) atoms into one Cu 
(group I) and one Tl (group III), forming ordered I-III- 
VI2 chalcopyrite compounds such as CuTlTe2 [H, US] ■ 

In Fig. 2(a) we plot the calculated band structure of 
HgTe under a e = 0.02 (001) tensile strain with Acf=70 
meV and CuTlTe2 in the chalcopyrite structure with 
Ac_F=76 meV. As we can see, a small gap is opened 
near the F point for both systems. Although the size 
is small, this anticrossing gap is protected by the lat- 
tice symmetry[21]. For the band structure calculation 
we employed density functional theory with a hybrid 
exchange-correlation functional, which can correctly pre- 
dict the band gaps of many zinc-blende and chalcopyrite 
semiconductors [22 - [24| . 

Comparing the band structure of HgTe under a e — 
0.02 (001) tensile strain and CuTlTe2, we find that the 
overall shape is very similar, especially near the band gap. 
In both systems, the s-like Fgc state falls below the p-like 
Tgy and Fy^ states, and the minimum gap occurs along 
the F — Xz line. This similarity between strained HgTe 



and CuTlTe2 indicates that the strain and chemical or- 
dering have the same effect in producing the crystal field 
splitting Acf at the top of valence band[26|. therefore, 
it could be an efficient way to tune the TI band gap. 

To achieve this goal, it is important to understand first 
how the splitting at the top of valence band is influenced 
by chemical ordering and what is the resulting depen- 
dence of the TI band gap. Based on the quasi-cubic 
modeljl^, [2^, and assuming the Fgc state is far away 
from the band edge, we know that the splitting of the 
Fg^ level into Fg^ and Fy^, under the tetragonal symmetry 
depends on two quantities: the spin-orbit splitting Ago 
and the crystal field splitting Acf. Acf is defined to be 
positive if the doubly-degenerate F51, is above the singly- 
degenerate F4„ state when the spin-orbit interaction is 
not considered, as shown in Fig. 1(a). For systems where 
Ago is much larger than Acf, the splitting between Fg, , 
and Fy^ is close to 2/3 of Acf- Previous studies |2]| 
on strained zinc-blende compound showed that the non- 
trivial gap depends on the sign and size of Acf: (i) when 
Acf < the gap occurs along of the F — Xx line near 
the F point and the gap increases quickly as a function 
of the magnitude of Acf; (h) when Acf > 0, the gap 
occurs along the F — Xz line near the F point and the gap 
increases slowly as a function of Acf- This can be seen 
clearly in Fig. 3(a), where the dependence of the non- 
trivial band gap on the size of Acf for HgTe is plotted. 
For example, when Acf=-100 meV, the gap is almost 40 
meV, but when Acf=100 meV, the gap is only 5 meV. 
The reason for the more significant gap increase with neg- 
ative Acf is that, F — Xx line has lower symmetry than 
F — Xz line, so the band anticrossing is more significant 
when the gap shifts to the F — Xx line. CuTlTe2 has a 
calculated Acf = 76 meV. This positive value explains 
why the gap shifts to the F — Xz line with only a small 
value of about 14 meV [Fig. 2(b)]. Based on this ob- 
servation, we know that large gap TI can only exist in 
zinc-blende derived compounds with large negative Acf- 

In Fig. 3(b) we plot the calculated Acf of CuAlTe2, 
CuGaTe2, CuInTe2 and CuTlTe2. As we can see, Acf 
increases from negative to positive as the group-Ill 
cations change from Al to Tl, i.e., from small light to large 
heavy elements. Considering that large negative Acf 
enlarges the non-trivial gap, one may intend to search 
compounds with small group-Ill cation as candidates 
for TIs. However, the requirement of band inversion at 
the F point excludes Al, Ga and In compounds because 
CuAlTe2, CuGaTe2 and CuInTe2 all have the normal 
band order. Their band gaps E^'P = E{rec) - S(Fg,y„) 
are all positive and decrease from Al to Ga to In com- 
pounds. This can be understood according to the band 
component of I-III-VI2 chalcopyrites: the Fgc state has s- 
like anti-bonding character localized on group-Ill cation 
and group- VI anion, whereas Fy^, and Fg^ states mainly 
have the p component of the group- VI anion hybridized 
with the d component of the group-I cationjl^, [2^. Two 
factors shift the level down from Al to Ga to In com- 
pounds [13, HO, H^: (i) the s orbital energy of Ga is 
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FIG. 2: (Color online) The calculated band structure along the high symmetry lines, Xx ; 2ZL(1 0) F : (0 0) — >■ Xz : 
^(0 1) of (a) HgTe with a (001) tensile strain and Acf = 70 meV, (b) CuTlTe2 and (c) Ag2HgPbSe4 at their equilibrium 
states. Xx and Xz are the notations of zinc-blende structure, and Xz corresponds to T in the chalcopyrite structure. Red 
and blue color are used to show the two spin-dependent bands clearly. 
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FIG. 3: (Color online) (a) The calculated non-trivial band 
gap as a function of Acf for HgTe. Here Acf is changed by 
tuning the (001) strain e. (b) The calculated Acf and (c) 
E'g-P of Cu-III-Te2 with III=A1, Ga, In, Tl. 



deeper than Al and (ii) In is much larger than Ga. For 
Tl, its s orbital energy, like Hg, is very deep due to the 
large relativistic effect, so its band gap is much lower 
than that of the corresponding In compounds. This is 
confirmed in Fig. 3(c), where we plot the calculated 
E'g-P of CuAlTea, CuGaTea, CuInTea and CuTlTca; only 
CuTlTe2 has negative Eg~P, i.e., its s-like Fgc state falls 
below the p-like states at F point (band inversion). Un- 
fortunately, CuTlTe2 has a positive Ac f = 76 meV, thus 
only a small Tl band gap of about 14 meV. 

In the above discussion, we have assumed that the Fgc 
state is deep inside the valance band, thus has no effect 
on the band splitting and the non-trivial gap of the TIs. 
However, if the Fgc is close to the band edge, then we have 
to consider its interaction with the band edge states. This 
is because when Aqf < 0, the band gap of the Tl at F 
point is between the unoccupied Fy^ and the occupied T^y 



(or Fgc, if it has a higher energy than Fgi,) derived state. 
The coupling between the T^y and Fgc states pushes the 
Tgy level up in energy, thus reduces the effective crystal 
field splitting between the Fy^ and Fg^ state and the non- 
trivial band gap. This is what we find for AgTlSe2 and 
AgTlTe2. According to our calculation, the non-trivial 
gap of AgTlSe2 is limited at F point with a very small 
size, 1 meV, although it has a large negative Acf = —50 
meV. Therefore, to reduce the interaction between the 
Tgy and Fgc states, one should move the Fgc level down, 
i.e., increase the magnitude of negative Eg~P as much as 
possible. 

The above analysis indicates that to design large 
gap chalcopyrite I-III-VI2 TIs, we face two contradic- 
tory requirements (large negative Acf and large neg- 
ative Eg~P). This severely limits the largest non-trivial 
gap obtainable for I-III-VI2 compounds. Through the 
direct calculation, we find that most of the already- 
synthesized I-III-VI2 have positive Eg~P and are normal 
semiconductors except CuTlSe2, CuTlTe2, AgTlSe2 
and AgTlTe2. But the non-trivial gaps of these four TIs 
are all small due to the positive Acf for CuTlSe2 and 
CuTlTe2, and small E'g-P for AgTlSes and AgTlTe2. 

To further increase the non-trivial band gap, we need 
to make both the Acf and Eg~P more negative. We 
find that this can be done by mutating two group-Ill 
cations in I-III-VI2 compounds to one group-II and one 
group-IV cation, thus forming the I2-II-IV-VI4 (I=Cu, 
Ag, II=Zn, Cd, Hg, IV==Si, Ge, Sn, Pb, VI=S, Se, Te) 
quaternary compounds. These compounds crystallize in 
either tetrahedral kesterite or stannite structures. Due 
to the increased chemical and structural freedom in the 
quaternary compounds, their band structure can be bet- 
ter tuned. Also, because they are structurally derived 
from chalcopyrites, their band structures maintain simi- 
lar characteristics [20!, ^295. Therefore, if these compounds 
have inverted band structure, they can also be TIs. 
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TABLE I: The calculated E^-p of l2-II-Pb-Vl4 (I=Cu, Ag, 
II=Cd, Hg, VI=S, Se, Te) in their ground-state structure. 
TM, TI and NI in the parentheses represent topological metal, 
topological insulator and normal insulator, respectively. 





Structure 


Te4 


Se4 


S4 


Cu2HgPb 


stannite 


-0.46 (TM) 


-0.32 (TI) 


0.07 (NI) 


Cu2CdPb 


stannite 


-0.21 (TM) 


-0.07 (TI) 


0.32 (NI) 


Ag2HgPb 


kesterite 


-0.37 (TM) 


-0.14 (TI) 


0.40 (NI) 


Ag2CdPb 


kesterite 


-0.12 (TI) 


0.18 (NI) 


0.72 (NI) 



Similar to the chalcopyrites, we need to have com- 
pounds that contain heavy group-IV elements so that the 
Fee level could fall below the Fg^ and Fy^ levels [20i [29|. 
Table I lists the calculated E'g-P of l2-n-Pb-Vl4 com- 
pounds. The results show that most of the Pb-Te and Pb- 
Se compounds have negative Eg~P at F and are therefore 
candidates for TIs. The calculation also shows all sul- 
phides and compounds containing other group-IV cations 
(Sn, Ge, Si) have positive Eg~P and are normal semicon- 
ductors. 

We first look at the band structure of Cu2HgPbTe4 
which has the most negative Eg~P. The overall shape 
near the F point is similar to those of ternary CuTlTe2 
and binary HgTe under (001) strain as shown in Fig. 2(a) 
and 2(b), indicating that the band structure character is 
kept in the cation mutation. However, Cu2HgPbTe4 is 
actually a topological metal (TM), because the conduc- 
tion band near L(N):^(0.5 0.5 0.5) point drops below 
VBM and crosses the Fermi level. The reason is that the 
conduction band state at L(N) point has similar charac- 
ter to the Fee state, so when the Fee energy is too low, 
the Lic{Nic) state energy is also below VBM, making 
the system metallic. To avoid this situation, therefore, 
we should search for TI material with mildly negative 
El-P. 

Our previous studyfH has shown that replacing 
Cu by Ag or replacing Te by Se can increase Eg~P, i.e., 
raising the Fgc and Lie energy level relative to Tgy and 
T-Yy, because (i) at the top valence band, the lower 4d 
level and larger size of Ag compared to Cu weakening 
the p-d hybridization, and the 4P level of Se is lower 
than 5p level of Te, which both shift the Fg,, and Fy^ 
levels down, (ii) the displacement of anion towards Pb 
in the Ag compounds and the smaller size of Se than Te 



also both reduce the Pb-anion bond lengths, increasing 
the energy of the Pb(s)-anion(s) antibonding states at 
the bottom conduction band. This expectation is sup- 
ported by the calculated band structure of Ag2HgPbSe4, 
which has no band crossing at the Fermi level and thus 
is a topological insulator, as shown in Fig. 2(c). Simi- 
larly, we predict that Cu2CdPbTe4 and Ag2HgPbTe4 are 
topological metals, while Cu2HgPbSe4, Cu2CdPbSe4 and 
Ag2CdPbTe4 are topological insulators. The results are 
shown in Table I. 

Among the four identified quaternary TIs, 
Ag2HgPbSe4 has the largest non-trivial gap of 47 
meV. This is because Ag2HgPbSe4 is more stable in 
the low symmetry kesterite structure with large group-I 
element, therefore, it has a large negative Acf (-51 
meV). It also has a reasonably large negative Eg~P gap, 
so the coupling between the Fgc and the Fg^, state is 
weak. Its band structure is shown in Fig. 2(c). As 
expected, we see the gap of Ag2HgPbSe4 occurs at a 
position along the low symmetry F — Xx line, consistent 
with our discussion above. 

In conclusion, we have shown that the non-trivial band 
gaps of zinc-blende derived topological insulators depend 
on the crystal field splitting at the top valence band as 
well as the size of the inverted s-p band gap. In general, 
a material with large TI band gap should have a large 
negative crystal field splitting and a moderate size of the 
inverted band gap. Compared to binary zinc-blende and 
ternary chalcopyrite compounds, these parameters can 
be more easily tuned through the chemical ordering in 
quaternary compounds. Based on this understanding, we 
have identified four ground state quaternary topological 
insulators, among which Ag2HgPbSe4 has the largest TI 
band gap of 47 meV because it has the optimal band 
structure parameters. 
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